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Abstract 

Scalar density cosmological perturbations, spectral indices and reheating in a chaotic inflationary 
universe model, in which a higher derivative term is added, are investigated. This term is supposed 
to play an important role in the early evolution of the Universe, specifically at times closer to the 
Planck era. 
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The most appealing cosmological scenario to date is the standard hot Big-Bang model. 
However, when this model is traced back to early times in the evolution of the Universe it 
presents some shortcomings ("puzzles"), such that the horizon, the flatness, the primordial 
density fluctuation problems, among others. 

Inflation [1] has been proposed as a good approach for solving the cosmological "puzzles" 
mentioned above. The essential feature of any inflationary model proposed so far is the 
rapid but finite period of expansion that the Universe underwent at very early times in its 
evolution. 

The present popularity of the inflationary scenario is entirely due to its ability to generate 
a spectrum of density perturbations which may lead to structure formation in the Universe. 
In this way, it has been mentioned that inflation is a good candidate for generating density 
fluctuations at scale significantly larger than the Hubble radius without breaking causality 
3]. In essence the conclusion that all the observations of microwave background anisotropies 
performed so far [^J support inflation, rests on the consistency of the anisotropies with an 
almost Harrison-Zel'dovich power spectrum predicted by most of the inflationary universe 
models. 

Among the different models we distinguish Linde's chaotic inflationary scenario j^J. This 
model allows various initial distribution of the inflaton field, providing a wide class of theories 
where inflation occurs under quite natural initial conditions. In this respect, it is possible to 
consider a model in which it starts at time closer to the Planck time, tp ~ 1CT 43 sec, where, 
it is supposed that quantum corrections to the effective Lagrangian play an important role. 
It has been suggested that the effective gravitational Lagrangian closed to the Planck era 
should include higher derivative terms expressed for instance, as a function of the scalar 

n 

curvature, R |5[. For simplicity, in this paper, we shall restrict to the case in which the 
effective action presents the i? 2 -term only. Generalization to include other more complicated 
terms looks straightforward. 

When the i? 2 -term of the effective action is varied with respect to the metric tensor 
g^v, it gives a term proportional to 1 H /1U in the field equations. This term coincides with 
that obtained as semiclassical limit of the interaction between a quantum matter field and 
gravity, via the expectation value of the corresponding energy-momentum tensor |6]. The 
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tensor 1 H flv is given by 

1 H fJiU = 2R. fl . l/ — 2g fll/ R. p ' p + 2RR flu — -g^ u R 2 . 

Because higher order curvature terms are suppressed by the inverse Planck mass squared, 
their effects will be significant near the Planck era. In this respect, a higher derivative 
term like i? 2 -term can be considered in addition to the Linde's chaotic inflationary universe. 
This kind of model as been treated in Refs. Q, Q. In Ref.Q the authors claim that the 
combined action of the R 2 term and the scalar field <f> may lead to Double Inflation i.e., two 
consecutive inflationary stages separated by a power-law expansion. However the analysis 
of perturbations made in j| is incomplete, because it was assumed that both contributions 



can be treated separately. These facts were elucidated in Ref.|8|, where double inflation 
was found to occur under the following conditions: the mass of the scalar particle has 
to be small compared to the scalaron mass m 2 <C M 2 , and the vacuum polarization has 
to dominates initially. Otherwise the combined action of both effects leads to a single 
inflationary phase. They also found a break in the spectrum of perturbations (the so called 
Broken Scale Invariant-BSI) for double inflation. 

The aim of the present paper is to study the effects that higher derivative term might 
have on inflation. Because (almost) all the literature on this topic works in the physical 
frame, leading to cumbersome expressions, we work here in the conformal frame where all 
the quantities can be easily derived. Here we study the single inflationary phase and left 
the double inflationary model for a future research |9j]. 

Specifically, we find that inflation can be enlarged by the action of the R 2 term. Actually, 
we concentrate our analysis in the case of a single inflationary phase driven by both contri- 
butions. We find, in agreement with previous work 3, that this inflationary phase occurs 
for a special set of initial values of the fields, constrained by the parameters of the model. 
We also determine analytically the scalar density perturbations produced during inflation 
valid for a non-separable effective potential, as the one obtained in this model, which has 
not been discussed in the literature. By using the observational constraint for the amplitude 
of density perturbations we find that m 2 /M 2 ~ 1CT 1 , showing consistency with a single 
inflationary phase. In Ref.ja] the authors use the ultrasynchroneous gauge allowing them 
to treat the problem in the physical frame, where it is found that longitudinal potentials 
are proportional to the ultrasyncroneous potential. Here we will consider longitudinal gauge 
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quantities, which are known to be gauge invariants Q| . Moreover, we obtain explicit analyt- 
ical expressions for the spectral indices when the universe leaves the horizon, showing small 
deviations from the Harrison-Zel'dovich spectrum, that can be observationally interesting. 
We also perform in this model the analysis of reheating, based on the modern theory of 
preheating. 

We consider the effective action given by 



S — d 4 Xy/^g 



(i) 



where R, as above, represents the scalar curvature, <fi is the "inflaton" scalar field, V(<f>) is the 
inflaton potential which drives inflation, and M is an arbitrary parameter with dimension 
of Planck mass, irtp ~ 10 19 GeV and k 2 = 87f/m 2 p . 

In order to calculate the scalar density perturbation we take the following conformal 
transformation 

= tt 2 (x)g^, (2) 

with 

fi 2 (x) = (1 + 3^), (3) 

where g^ u and g^ represent the original and the transformed metric, respectively. Introduc- 
ing a new field ip defined by 

^h Hl + ip ] - (4) 



where j3 = a/2/3, and substituting eqs. (j2EJ) in the action (0), it becomes a Hilbert- 
Einstein-type action given by 

s = J dv=i{^ + ~w + \e-^m 2 + 

(5) 

where the potential V(4>, ip) takes the form 



V(M) = 7V(l-e"^) 2 + e~ 2 ^V(<p). (6) 

4 K z 

Notice that, if we disregard the inflaton field 0, then inflation is described entirely by gravity, 
since the model becomes supported by the i? 2 -theory is] ]. 

Let us introduce into the action (jHJ) the conformal Friedmann-Robertson- Walker (FRW) 
metric 

ds 2 = dt 2 - d(i) 2 dxidx\ (7) 



where dxidx 1 and a(i) represent the flat three-surface and the conformal scale factor, re- 
spectively. Considering also that the corresponding scalar fields <fi and ip are homogeneous, 
i.e., = <j)(t) and ip(t), then, after variations of this action respect to the different field 
variables, the following Equations of motion are found 



3H 



K 



2 V 2 



2 + V{M) 



ip + 3Hip 



(3k 



and 



+ 3H<p = (3nip(p- e 



dip 



(8) 



where the dots denote differentiation respect to the time i and H = a/a defines the Hubble 
parameter in the conformal frame. 

In the slow-roll over approximation the field eqs. (|SJ) become, 



3H' 



3HiP = - 



dip 



(9) 



and 



3Ht 



M dV{iPA) 



This set of eqs. is valid if the following conditions are met 



max{ip 2 , < V(ip, 



and 



(V^,0), ) 2 « 



-^V{iPA),H 



where V (ip , <p) ^ and V(ip,<p) t ^ represent the derivatives of the potential V(ip,(p) respect to 



the fields <p> and ip, respectively. A so. 



inflationary potential, V((p) = \m 2 <p 2 {% and M 2 ^ 2m 2 is given by 



ution to the equation of motion (JHJ) for the chaotic 



where, we have taken the integration constant C = =>- e /3K ^° — 1 = ^^-^n 2 ) The 
subscript denotes the values of each field at the beginning of inflation. The approximation 
involve in (fTUJ) is used for simplicity. In this case, we find 



(M 2 - 2m 2 ) 
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H K 2 



and 



where 



if) = —In 



SH 



(11) 



(12) 



- 2 _ M 2 _ -/3ki/>o \2 _j_ rn 2 K 2 (f) 2 2/3^0 _ k2 Vo 

1 (i "~ 3 ' 

Since we are mainly interested in the existence of inflationary stage, we should take the 
case in which M > y/2m. Here, the shape of the potential exhibits an appropriated form 
for inflation to occur. This is not the case if M < V2m, since the potential presents an 
unacceptable peak for certain values of the scalar fields, which would force us to restrict 
strongly the range of these fields for inflationary models to occur. 

The period of the superluminar expansion can be determined by the value of the inflaton 
field where the slow-rollover approximation breaks down. This happens when |0| ~ \3H<p\. 
This condition gives 

<2,-^ j2 (M 2 -2m 2 ) 

♦"^'-W < 13 > 



o 



and the inflationary period results in 

t ~ L (14 ) 

e " (M 2 -2m 2 ) QH ' { ' 

where the subscript e denotes the values at the end of inflation. 

The requirement for solving some of the puzzles present in standard Big-Bang cosmology 

is 

a(0) 5(0) ~ 1 } 
This results in a condition for the initial value of the inflaton field given by 

9 (3 2 H 2 eg k 2 > 8 2 [M 2 - 2m 2 ] (16) 



or 

M 2 1 



O r^u 



M 2 



46 - 



(3 2 m 2 



(17) 
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which dependents on the parameters M and m, as we could expect for two-field models with 
a no separable potential. 

In order to calculate the scalar density fluctuations, we follow a similar approach to that 
described in ref. |l4j. We take the perturbed metric in the longitudinal gauge, which is 
known to be manifestly gauge- invariant Q|. Therefore, in comoving coordinates we write 

ds 2 = (1 + 2$)dP - a 2 (l - 2$) 8 i:i dx l dx j , (18) 

where the scalar metric perturbation fields are functions of the space and time coordinates, 
i.e., <J = and ^ = \i>(f,t). 

In order to describe a linear theory for the cosmological perturbations, we introduce 
fluctuations into the scalar fields <fi and ip, so that we write <j)(x, t) = <p(t) + 5(p(x, t) and 
ip{x,t) = ip(t) + 5ip(x,t), where the background fields, (fi(t) and i[>(t), are solutions to the 
homogeneous Einstein equations, and the deltas are small perturbations that represent small 
fluctuations of the corresponding scalar fields. 

Since the spatial part of the resultant energy-momentum tensor becomes diagonal, then 
it found that $ = 1 l r . With this result, we can write the perturbed Einstein field equations, 
in which we can take spatial Fourier transform for the different variables. It is found that 
each mode satisfies the corresponding linearized perturbed field equation, given by 

k 2 r • -i 

' ~^cf>S(f) + #V , (19) 



<5> + H<5> 

2 



and 



2(0 + 3i/0)l> + 40<l -2/?k0^$, (20) 



/ 2 o2 2 

6$ + 3H5xP + (- - + V^W 



+(3Ke- p ^<p5<p + V^5<j) 



2(^ + 3^)$ + 4^$ + /?Ke- /3K ^0 2 $. (21) 
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If in this set of equations we neglect terms proportional to $ and two time derivatives, and 
we consider non-decreasing adiabatic and isocurvature modes on large scale, i.e. k <^aH, 
the following equations are obtained j3| 

.2 



$ = -= 

2H 



(22) 



3H6<P + (e^V^)^ + e^V^Scj) + 2e^V^ = 0, 



and 



(23) 
(24) 



where V is the potential specified by eq. ©• 

For V((f>) = m 2 2 /2 we can solve the above equations, following a similar approach to 
ref. 1 15]. It is found that 



and 



where 



and 



1 



-n 



— = i[^i + c 2 / 2 (^,0)] 



V 2 V 2 

v 2 + v 2 1 



/l 



a 



(2 



3M 2 /o m 2 2 
4k 2 



■(2-e 



2C e i 



(25) 



(26) 



(27) 



h = h(lf>,4>) = a 2 (0 e ,Ve) + 



3M 2 



X 



m 



C e2 I — c^e 



2-PKip 



V In 



_(3M 2 (e"W - 1) + 2K 2 m 2 2 e-^)_ 

+ C e i(e-^-l) 



3M 2 (e" /3 ^ - 1) m 2 2 e~^ 



4k 2 



where «i and a 2 are arbitrary constants and G e \ and C e2 are given by C e \ 



v. 



__(l_ c ^ ) + 



2^2 



and C e2 = respectively. 



In the expressions f|25H27|) . terms in proportion to C\ and C 3 represent adiabatic and 
isocurvature modes, respectively. 
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^From eqs. (|25j) and (|26j) the curvature perturbation due to primordially adiabatic 
fluctuation at the end of inflation, when the scale crosses the horizon is given by 

8_p 
P 



^[1 -/(**)] + 



(28) 



where 



h 



In the limit M~ 







[fx + h ~ IJ 

■0 = 0, which gives / = Cte. and thus we obtain the standard re- 
sult for the density perturbation Sp/p oc H 2 /4>- The magnitude of the density perturbations 
take the form 

8_p 
P 



[M 2 - 2m 2 ] 



x 



1 + 



2m 



. , . (29) 

v/(M 2 - 2m 2 ) f 

where <ph represent the field scalar at witch the universe crosses the Hubble radius and 
f t f >h is the function / evaluate at fa. Since a(th)/a(t e ) = kH~ x we find 



(M 2 - 2m 2 ) 
H 2 ^ 2 



HHok- 1 ) + 



6 



(30) 



and in view that the log term is of the order of O(10 2 ) and 5p/p ~ 10" 5 , we find, from 
eqs. f)29p and (|3(jp . the following order for the parameters of our model m 2 /M 2 ~ 10 _1 . 
For this order the initial value of <fi becomes restricted from below n<p > 5. 

Another interesting quantity to investigate is the spectral index ns of the power spectrum 
of primordial fluctuations, which becomes specified by P(k) oc k ns . It is well known that 
most of the inflationary scenarios are compatible with an almost Harrison-Zel'dovich spec- 
trum, i.e., ns ~ 1, and thus, it provides detailed information about the early time in the 
evolution of the universe. This prediction of the inflationary universe models is sustained 
by the COBE data [3j], which are statistically compatible with the ns ~ 1 result. 

In our model, we could obtain the spectral index by using the following standard expres- 
sion 

din 



n s - 1 



5_p |2 
P 



d In k 



(31) 



where, together with eqs. (|2T?|) and (|3Ti|). it is obtained, 



-2 



[\n(H 2 k-i) + i 



(/* - 1) + 



(M 2 -2m 2 ) 



fk 



fk) 2 + 



(M 2 -2m 2 )Jk 



where fk represents the function / expressed in terms of k and f' k = dfk/dlnk. The right- 
hand-size of the latter equation represents the deviation from the fig — 1 Harrison-Zel'dovich 
spectrum. For instant, if we use the constraint obtained aboved, i.e., m 2 /M 2 ~ 10 



-i 



we 



obtaine that fig ~ 0.98, which differs from the Harrison-Zel'dovich (HZ) value approximately 
in a 2%. We should note that this deviation is inside of the observational range, since it 
has been reported that the scalar perturbations do not differ from scale invariant by a large 
amount, since \ns — 1| ~ (9(0.1), i.e., inside of the 10% j^J. 

In addition to the scalar curvature perturbations there are tensor or gravitational wave 
perturbations. These perturbations can also be generated from quantum fluctuations during 
inflation. Since we have chosen to work in the Einstein conformal frame we can use the 
standard results for the evolution of tensor perturbations. The spectrum of the gravitational 
wave Pq is known to be g iven by Pq = 8k 2 (H/2it) 2 , and where the spectral index is defined 
by uq = dlnPc/dlnk |l0l 121]. Thus, we find in our case 



n G 



3M 4 



ln(H k 



\n(H k 



+ 6 



where (3 = 9M 2 (M 2 — 2m 2 ) / (8/€ 4 m 2 ). We could give an order of magnitude of this parameter 
if we use the value obtained for the quotient m/M determined above. Our result gives 
fiQ pa —0.22. Notice that this value differs in about a 20% of that value predicted by the 
HZ spectrum which gives tiq ~ 0. Thus, our model presents a significant deviation from the 
HZ result. Only astronomical observations will tell us which spectrum will be the right one. 

Note that every expression was determined in the Einstein conformal frame. In order to 
see what happens with the same expressions but in the physical frame, we have to use the 
original conformal transformation eqs. (J2J) and (J3J). The quantities in the Einstein frame are 
related to the physical frame by means of a conformal transformation given by eq.(J2J). In 
both frame, the gauge-invariant perturbations are related by |2J 



- SQ 
$ 



(32) 



where from eqs. (3) and (4) it is found that 5Q = ^Q8ip/2, with 5ip 
In the physical frame, the gravitational potential becomes 

1 



H. 



$ = C x 1 



ail 2 



ail 2 dt' 



(33) 



Comparing both expressions, it is observed that they are very different, but, during inflation, 
where the Hubble parameter remains practically constant, we expect that the quadratic 
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scalar curvature term in the effective Lagrangian varies slowly, so that the f2 conformal 
factor could be considered constant during this period. Therefore, the adiabatic fluctuations 
become described by the same expression in both frames. The same argument can be draw 
for the spectral indices n$ and since it depends directly of the expression for 5p/p and 
H respectively. 

Every inflationary model needs a mechanism to reheat the universe, i.e., a physical mech- 
anism capable to patch the inflationary era with the radiation dominated phase (needed 
for example in nucleosynthesis). The current understanding of this process comprises three 
phases: preheating (the resonant amplification phase), the perturbative decay and thermal- 



ization 



11 1 . Reheating in the R 2 model has been considered by Suen and Anderson |23[ by 



24( have studied 



solving the semiclassical backreaction equations. Recently Tsujikawa et al. 
the preheating phase with nonminimally coupled scalar fields. However, as far as we know 
the study of preheating in the model that we are considering here, has not been discussed 
in the literature. 

In order to study the preheating phase we consider the inflaton field coupled with another 
scalar field, say x through the interaction term 

-^W, (34) 

where g is the coupling constant. During the coherent oscillations of the inflaton field around 
the minimum at the origin, which is precisely the period where reheating occur, the scalaron 
field moves towards the minimum, but in a time scale shorter than that of the inflaton </>. 
In essence we can consider it nearly constant during preheating; which is the first part of 
reheating. Under these assumptions the inflaton field equation reads 

4> + 3H4> = -e^(m 2 + g 2 ( X 2 ))4>, (35) 

where we have replaced x 2 D Y the quantum expectation value, (x 2 ), doing manifest the 
mean field approximation involved. The precise definition of this object will be given in 
short. Furthermore, the equation of the x field is expressed by 

X + 3H X - V 2 X = -e^V0 2 )x- (36) 



As we can notice here, the preheating phase proceed very similar to the standard way, 
because the only difference is the exponential damping term in the coupling. By considering 

11 



the field x as a quantum operator we can expand it as 

d 3 k 



X 



a kX k(t)e~^ + a{xl(i)^ , (37) 



(2vr)3/2 

where the a's are annihilation and creation operators, and the mode functions Xk(t) satisfy 
the equation 



Xk + 3Hxk + 



k^ 2 



+ e^g 2 cj) 2 



Xk = 0. (38) 



The expectation value of \ is expressed in terms of the mode functions as 

(X 2 ) = A I dkk 2 \x\\. (39) 



2vr 2 j 

Eq. ()38|) can be rewritten as a harmonic oscillator with a time dependent frequency 

x k + n 2 k (t)x k = o, (40) 

where we have defined the variable X k = a(t) 2 ^ 3 Xk and the frequency is given by 

n J=(i)' + «*W < 41 > 

Usually the last two terms in (|41|) are of the same order and also very small compare with 
the effective mass m x , so we can neglect these terms as a first approximation. In this case 
the frequency depends on the scale factor d(t), the scalar field <f>(t) and if>(t). 

During the oscillations of the inflaton field 0, the scalaron if) first rolls towards the mini- 
mum and after that start to oscillate around it, right after a few oscillations of </>. Preheating 
occurs here before if) start to oscillate. If we define ip Q as the value at this point, an approx- 
imated solution of (jSSJ) is 

(p(t) ~ $(t)sin(wt), (42) 

where to = e^m and the amplitude $(f) changes slowly with t. Inserting this solution in 
eq. ()41|) we find a Mathieu type equation. It controls both the narrow and broad resonance 
particle production. We have performed numerical studies of that equation finding no sur- 
prises, just the evident effect of damping terms in (|35|) and (J4(Jj) . After if) start to oscillate 
(for if) < if) ), the system becomes completely dominated by if), implying another period of 
reheating, this time without any resonance amplification, in which the if) field decays. A 
detailed research of preheating in these models, considering both inflationary realizations, 
is under study 0. 
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